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Abstract 



We calculate perturbative renormalization properties of the topological 
charge, using the standard lattice discretization given by a product of twisted 
plaquettes. We use the overlap and clover action for fermions, and the 
Symanzik improved gluon action for 4- and 6-link loops. 

We compute the multiplicative renormalization of the topological charge 
density to one loop; this involves only the gluon part of the action. The 
power divergent additive renormalization of the topological susceptibility is 
calculated to 3 loops. 

Our work serves also as a test case of the techniques and limitations of 
lattice perturbation theory, it being the first 3-loop computation in the 
literature involving overlap fermions. 

Keywords: Lattice QCD, Topology, Lattice perturbation theory, Over- 
lap action, Improved actions. 
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Introduction 



Topological properties of QCD are among those most widely studied on the lattice. 
Various methods have been used to this end, involving renormalization, cooling, fermionic 
zero modes, geometric definitions, etc. In recent years, the advent of fermionic actions, such 
as the overlap, which do not violate chirality, has brought a new thrust to the subject. 

In this work we compute the renormalization constants which are necessary in order to 
extract topological properties, in the "field theoretic" approach, from Monte Carlo simula- 
tions using Wilson or Symanzik improved gluons, and clover or overlap fermions.In 
particular, we compute the multiplicative renormalization Zq of the topological charge den- 
sity, to 1 loop in perturbation theory and the power divergent additive renormalization 
M(g 2 ) of the topological susceptibility, to 3 loops. 

The main motivations for doing this work are: 

• To enable comparison between different approaches used in studying topology, so that 
a coherent picture of topology in QCD may emerge. 

• To enable studies, in numerical simulations, of quantities involving the density of 
topological charge, q(x), rather than only the integrated charge; this is necessary, e.g., 
for studying the spin content of nucleons. 

• As a feasibility study in lattice perturbation theory: Indeed, this is the first 3-loop 
calculation involving overlap fermions. 



Computation of Zq 

Our first task is to compute the the multiplicative renormalization Zq [1] of the topo- 
logical charge density q L (x) to one loop, using the background field method. We use the 
standard definition of q L , given by a product of twisted plaquettes 

1 ±4 

Ql( x ) = —w-2 E e^ pa TT[U^(x)U pa (x)} (1) 

Z 71 t ivpa=±l 

{t-n,v,p,o = —£n,v,p,(r)- n Miy (x) is the parallel transport matrix along a 1 x 1 Wilson loop; in 
standard notation 

IWx) = U,(x) U u (x + n) Ul(x + u) Ul(x) (2) 

The classical limit of the operator shown in Eq.(l) must be corrected by including a 
renormalization function Zq, which can be expressed perturbatively as 

Z Q = 1 + Zi • g 2 + ■ ■ ■ , Zi = Z n • N c + Z 12 • -J- (3) 

We perform a calculation of Z\\ this involves only the gluon part of the action. 
In the background field method, link variables are decomposed as 
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U,(x) = V,(x)U c ,(x) (4) 
in terms of links for a quantum field and a classical background field, respectively 

V^x) = e l9Q ^ x \ U CfM (x) = e mB ^ x) (5) 
The N c x N c Hermitian matrices Q M and can be expressed as 

Q fl (x)=t a Q^(x), B II (x)=t a B;(x), Tr[t a t b ] = ±5 ab (6) 

The perturbative nature of our calculation requires a choice of gauge. An appropriate 
gauge-fixing term is 

S 9 /=^EE Tr\D-Q»D~Q v \ (7) 

This term breaks gauge invariance with respect to Q^, as it should, but succeeds in keep- 
ing the path integral as a gauge invariant functional of B^. The definition of the lattice 
derivative, which is covariant with respect to background gauge transformations, is 

D-{U c )Q v (x) = U^(x - e^)Q v {x - ejU^x - e„) - Q u (x) (8) 

The diagrams involved in the one- loop calculation of Zq are shown in Figure 1. 




Fig. 1: Diagrams contributing to Z\. The black bullet stands for topological charge vertices. 

External lines correspond to background fields. 



We use the Symanzik improved gauge field action, involving Wilson loops with 4 and 6 
links. In standard notation, it reads 

2 r 

C ReTr (1 - f/pi aqU ette) 

plaquette 

+ Ci ReTr (1 - f/ rcc tanglc) 

rectangle 



+ c 2 ]T ReTr(l-f/ chair ) 



chair 



+ C 3 ^ Re Tr (1 — f/p ara llclogram ) 
parallelogram 



(9) 
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The lowest order expansion of this action, leading to the gluon propagator, is (see, e.g., 
Ref. [2]) 



2 J-Tr/a (2tt) 4 ^ ^ ; 



(10) 



where: 
and: 



1 - V) [Co - Ci a 2 ^ 2 - C 2 a\kl + fc*)] 
2 . ak h 



k 



— sin ■ 



'/i 



^ 2 = E^. 



a 2 

The coefficients Cj are related to the Symanzik coefficients q by 

Co = c + 8ci + 16c 2 + 8c 3 , Ci = c 2 + c 3 , C 2 = c x - c 2 - c 3 

The Symanzik coefficients must satisfy: c + 8ci + 16c 2 + 8c 3 = 1, in order to reach the 
correct classical continuum limit. 

Our calculations are performed without any assumptions on the values of the external 
momenta p±, p 2 : This is safest for the topological charge operator, otherwise one may easily 
end up with indeterminate expressions. 

For our purpose, we must express all potentially divergent integrals / in terms of a con- 
tinuum counterpart I cont (evaluated in D = 4 — 2e dimensions), plus all lattice contributions 
Ibatt- The latter are the ones which will determine Zq. The three expressions shown below 
form a basis set for all the divergent integrals encountered in this calculation; we point out 
the need to handle also a 3-point form factor [Cp V below) [3] 



B(a,p) 
Bp(a,p) 
C ilv {a,p 1 ,p 2 ) 



(ka) 2 " f d D k 1 

J (2*)° V (k + ap)* ~ {P) + ° Lat 
(ka) 2e f d D k sink,, 



Ka) c r a k 



(2tt) d k 2 (k + ap) 2 
(ka) 2e f d D k sinkp sink v 



Bp(p) + h Latt 



(11) 
(12) 



a J (2tt) d k 2 {k + ap 1 ) 2 {k + ap 1 + ap 2 ) 2 



Cp V {pi,p 2 ) + 



QAtt 2 v 



— In n 2 a 2 — In An 



+(Air) 2 P 2 + ~f E -2ir 2 P 1 



(13) 



where is Euler's constant, B, Bp, Cp U are the continuum counterparts of B, Bp, Cp U , 
respectively; Pi, P 2 are [4] 

Pi = 0.15493339023109021(1), P 2 = 0.02401318111946489(1) 
and: In r ^— rr^ f - - - In An - Inn 2 a 2 + + P 2 (14) 
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(4tt 


f( 


e 


^Latt 


= pA 


Pi 
16 


1 7 
~2 V 



As we see, Iiatt contains poles in e; indeed, diagrams (c) and (d) of Figure 1, taken separately, 
exhibit such poles ( (d) oc — 1/e — lnn 2 a 2 ). These cancel, however, upon summation, as is 
expected by the fact that Q does not renormalize in the continuum. 

The complicated algebra of perturbation theory was carried out using our package in 
Mathematica. The calculation of Zq is particularly involved in the present case, involving 
propagators and vertices from the improved gluonic action. In particular, the calculation 
of diagram (d) involves a summation of > 1 000 000 different algebraic expressions at 
intermediate stages. 

Our results for Zq are listed in Table I. In all calculations that involve the parameters 
q, we choose a standard set of values as in Ref. [2]. The choice of the sets of parameters 
correspond to the most popular actions: The first set corresponds to the plaquette action, 
the second set corresponds to the tree-level Symanzik improved action [5] and the next 6 
sets correspond to the tadpole improved Liischer-Weisz (TILW) action [6,7] for 6 values of 
beta 

P = 8.60, 8.45, 8.30, 8.20, 8.10, 8.00 

The last two sets correspond to the Iwasaki [8] and DBW2 [9] actions respectively. 

In the case of the plaquette action, our result agrees with the known result of Ref. [1]. 

It is worth noting that the value of Z\ (and of e 3 in Table II) for the DBW2 action is the 
smallest one, leading to a renormalization factor Zq closer to 1 (and M(g 2 ) closer to 0, see 
below. This would single out the DBW2 action as a better candidate for studies 
of topology. 

Computation of M(g 2 ) 

The second task we attend to is the calculation of the power additive renormalization of 
the topological charge susceptibility, which is defined as 

X L = E<<7,(^(°)> ( 15 ) 

X 

X L develops an unphysical background term which becomes dominant in the continuum limit 

Xl (9 2 ) = ^Z Q {g 2 )\ + M(g 2 ) (16) 

The power divergent additive renormalization of x L can De written perturbatively as 

M(g 2 ) = e 3 -g 6 + e 4 -g 8 + --- (17) 

We first compute the 2-loop coefficient e 3 . Its dependence on the number of colors has the 
form: 

e 3 = (iV c 2 -l)iV c e3,o (18) 

This result is evaluated for several sets of values of the Symanzik improvement coefficients. 
Figure 2 shows the diagram contributing to . 
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Fig. 2: Two-loop diagram contributing to e%. Black bullets stand for topological charge vertices. 

The 3-loop term e 4 of the expansion of M(g 2 ) equals 

e 4 = el + e{ (19) 

where e{ stands for the fermionic contribution to e 4 (csw is the coefficient in the clover 
action) 

e{ = N f (N* - 1)N C • (e 4)0 + e 4 ,i c sw + e 4 , 2 4 W ) (20) 
and el is the purely gluonic contribution, which is expressed as in Ref. [10] 

e 9 4 = ^(N 2 C - l)(1.735iV c 2 - 10.82 + 73.83/iV c 2 ) x 10~ 7 (21) 

In fact, what we are interested in, is the calculation of the parameters e 4t o, e^i, e^. This 
task is performed using both overlap and clover fermions (Clearly, overlap fermions 
involve only the parameter e 4i0 ) • Figure 3 shows the 3-loop diagrams contributing to the 
evaluation of e{. 





(a) 

Fig. 3: Diagrams contributing to e A . Straight lines correspond to 
overlap or clover fermions and wavy lines correspond to Wilson gluons. 

The propagator and vertices for overlap fermions can be obtained from the following 
expression of the overlap action 

^Overlap O 

Y, V(n)D N (n,m)V(m) (22) 

n, m 

where D N is the massless Neuberger-Dirac operator [11] 

M is a real parameter corresponding to a negative mass term. M must lie in the range 
< M < 2r, r being the Wilson parameter (in our case r = 1). X is a Hermitian operator 
which can be expressed, in momentum space, in terms of the Wilson-Dirac operator Dw 
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X(q,p) = a [ D w (q,p) - — 



= aX (p)(27r)^(q -p) + aX l {q,p) + aX 2 (q,p) + 0{g 6 ) 



(24) 



where : A (p) = - 7 M s m a ^ + - zJ 1 ~~ cos a W 



X 1 (q,p)=gJ d 4 k5(q-p-k)A,(k)V 1 ^p + ^) 
with : Vi^(g) = 27^ cos ag M + r sin aq^ 

and: X 2 (q,p) = y ^ - J^^ Kq ~ V ~ h - k 2 ) 

A,(k 1 )A,(k 2 )V 2 Ap+Y + Y ) 
with : V2 iAt (g) = —i^^a sin ag M + ra cos ag^ 

The clover (SW) fermionic action [12], in standard notation, reads 

/ 



^ = ^ E T r [1 - H^(x)] + £ £(4r + m)^f(x)f/; f (x) 
y x,n,v fx 



1 



f X,fl 



+ip f {x + fi){r + 7/1 ) U^ix^fix) 



(25) 
(26) 

(27) 



(28) 



where : F M „ = ^ (Q^ - Q^) (29) 
and . Qfiu U Xj x+^Ux-^n^ x+n+vU x +n-\-v, x+vU x +v, x 

~\~ U Xj x+i/L ^x+v, x+v—ijJJx+v—h, x—fjJJx—^i, X 
~\~ U x , x-^Ux-fj,, x—ix—uUx—n—v, x—vUx—v, X 

U x ,x—vUx—v,x—v+ijJJx—v+ix,x+ijJJx+ix,x (30) 

The clover coefficient Csw is treated here as a free parameter; r is the Wilson 

parameter; / is a flavor index; = (i/2)[y fM , 7J. Powers of the lattice spacing a have been 
omitted and may be directly reinserted by dimensional counting. 

In performing this calculation, a large effort was devoted to the creation of an efficient 
3-loop "integrator", that is, a metacode for converting lengthy 3-loop integrands into 
efficient code for numerical integration. Some key features of the integrator are listed in the 
Appendix. 

Table II contains our results for e 3 (cf. Eqs.(17, 18)) for different gluonic actions. 
Fermions do not contribute here. Our results for the case of the plaquette action agree 



7 



with older known results (see, e.g., [13,14]). Tables III and IV list our results for e 4 , using 
the clover and overlap actions, respectively. 

Figure 4 shows the coefficients e 4 , e 4j i, e 4j2 of the clover result for different values of 
the bare fermion mass m. For ease of reference, Figure 5 presents the total values of e 4 for 
different choices of csw, with Nf — 2, N c = 3. 

Figure 6 exhibits the dependence of e 4j o , using the overlap action, on the parameter M . 
The total value of e 4 in this case is shown in Figure 7. 



Discussion and Conclusions 

We have calculated in this work the perturbative renormalization coefficients for the 
topological charge and susceptibility. The topological charge, defined via a usual "twisted" 
product of plaquettes, was renormalized multiplicatively to one loop using a variety of gluon 
actions presently employed in numerical simulations. The details of the fermionic action 
only enter this computation at 2-loop order. The 1-loop coefficient is rather pronounced for 
most actions but it is quite suppressed for the particular case of the DBW2 action. 

For the topological susceptibility, we performed a calculation of its power-divergent, addi- 
tive renormalization M(g 2 ) to 2 and 3 loops. Fermions do not contribute to the 2-loop result. 
Among the various gluonic actions, the DBW2 action leads again to a coefficient which is 
quite suppressed, by more than 2 orders of magnitude as compared to the standard Wilson 
action. Based on these results, the DBW2 action is singled out as the favourite candidate 
for simulations involving measurements of topology in the "field theoretic" approach. 

The 3-loop calculation was performed in the presence of clover, as well as overlap 
fermions, and Wilson gluons. Indeed, this is the first 3-loop calculation in the literature 
involving overlap fermions. To facilitate a comparison of various contributions to M(g 2 ), we 
write them in the particular case of N c = 3 with Wilson gluons: 

M(g 2 ) = g 6 1.655 x 1(T 5 + g 8 6.50 x 1(T 7 + 24 g 8 N f x (31) 

The first two summands above are clearly the purely gluonic contributions, at 2 and 3 
loops. The value of x can be read from Table III for clover fermions (x = e 4j0 + e 4 i csw + 
e 4 ,2 c| w ) and from Table IV for overlap fermions (x = e 4i o). In all cases, fermionic contri- 
butions are of the same order of magnitude as 3-loop gluon contributions, but still only a 
small fraction of the 2-loop result. 

Our results can be used to enhance a number of related computations. In some cases, 
the local definition of q L (x) is used (e.g., [15]), with renormalization estimates coming from 
simulations; other investigations propose non-ultralocal definitions of q(x) (e.g., [16]), which 
clearly are more expensive to simulate, but can circumvent renormalization; in yet other 
cases, operational/numerical definitions of q(x) are proposed (e.g., [17]). Even for the inte- 
grated topological charge Q, where fermions with exact chiral symmetry offer us in principle 
unambiguous ways of extracting Q from lattice configurations (see, e.g., Refs. [18,19]), a 
comparison with other time-tested definitions is called for. In all cases, it would be very 
important to verify that a consistent picture of topology in QCD emerges from the various 
approaches. 
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Appendix 



Here we briefly describe some features of our "integrator" program; this is metacode, 
written in Mathematica, for converting lengthy 3-loop integrands into efficient Fortran code 
for numerical integration on finite lattices of different sizes L 4 . The results r(L) are subse- 
quently extrapolated to L — > oo, via a fit to a large class of functions containing powers of 
L and InL, as is a priori expected; systematic errors on r(oo) are also produced as a result 
of these fits. 

A number of optimizations are implemented, as listed below; the resulting output code 
runs many orders of magnitude faster than output of automatic code conversion programs, 
such as, e.g., the built-in converter of Mathematica. 

The input expression to the integrator is an integrand which depends on 3 four-momenta 
Pi > Vii Pz ■ It has the form of a sum of terms (typically tens of thousands); each term is a 
product of trigonometric functions of combinations of Pi , possibly including also propagators 
for Symanzik improved gluons and overlap fermions. 

The following optimizations take place: 
A number of time consuming common ingredients are precalculated numerically and stored, 
for all possible momentum values on a finite lattice; such ingredients are trigonometric func- 
tions, the overlap propagator, and the Symanzik propagator. While an expression in closed 
form exists for the latter [2], in practice it is considerably faster to invert the propagator 
numerically and store it; after all, in either case the final results cannot be presented as 
analytic functions of the parameters q , since the dependence on q is not polynomial (unlike 
the case of csw)- 

The symmetries of the integrand are exploited to reduce the volume of the integration region, 
with due attention paid to correct counting of points at the borders. 

All 3-loop diagrams, with the exception of those in the form of the Mercedes emblem, 
contain two non-overlapping loops (i.e., loops with no shared propagators); consequently, 
their integrands can always be written (after some trigonometry) as a sum over expressions 
of the form: 

/ dptfM) (/ dp 4 2 f 2 ( Pl ,P2)) (J dptf 3 ( Pl ,p 3 )) (32) 

Integrations over p 2 and p% can then be performed sequentially, rather than in a nested fash- 
ion (actually, we perform them simultaneously, see next item), thus reducing their execution 
time to that of a 2-loop integral. 

The integrand is organized as an inverse tree, i.e. summands having the same functional 
dependence on the innermost integration variable are grouped together, and the innermost 
integral is performed once for each group; the procedure is then iterated for outer integrals. 
In practice, this procedure can save one or two orders of magnitude in execution time of the 
innermost integrals, which are the most expensive ones. 

Terms with different polynomial dependence on unspecified parameters such as N c or csw, 
etc., but with otherwise similar functional form are treated simultaneously, and the result 
is presented as a polynomial in these parameters. Also, for parameters whose values are 
read on input (masses, Symanzik coefficients), the code runs in parallel for different sets of 
values, in order to avoid computing the same quantities several times. 
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TABLES 

TABLE I. The values of Zn and Z12 (Eq.(3), Figure 1), for various values of the Symanzik 
coefficients cq, ci, C2, C3 (02 = 0). 



Action 




CO 




Cl 




C3 


Zu 




Z\2 


Plaquette 




1.0 














-0.33059398205(2) 





2500000000(1) 


Symanzik 




1.6666666 


-0 


083333 








-0.2512236240(1) 





183131339233(1; 


TILW, & = 


S.60 


2.3168064 


-0 


151791 


-0 


0128098 


-0.20828371039(3) 





147519438874(3; 


TILW, = 


S.45 


2.3460240 


-0 


154846 


-0 


0134070 


-0.20674100461(1) 





146259768983(1; 


TILW, (3 = 


S.30 


2.3869776 


-0 


159128 


-0 


0142442 


-0.20462181183(1) 





144531861677(4; 


TILW, /3 = 


S.20 


2.4127840 


-0 


161827 


-0 


0147710 


-0.20331145580(1) 





143464931830(1; 


TILW, = 


S.10 


2.4465400 


-0 


165353 


-0 


0154645 


-0.20162651307(1) 





142094444611(2; 


TILW, p = 


S.00 


2.4891712 


-0 


169805 


-0 


0163414 


-0.19954339172(1) 





140402610424(1; 


Iwasaki 




3.648 


-0 


331 








-0.15392854668(1) 





105132852383(2; 


DBW2 




12.2688 


-1 


4086 








-0.0617777059(4) 





038277296152(6; 



TABLE II. Evaluation of e 3i0 (cf. Eqs.(17, 18), Figure 2) with Symanzik im- 
proved gluons, for various values of the coefficients cq, ci, C2, C3 (02 = 0). 



Action 




CO 


Cl 




C3 


e 3 , x 10 7 


Plaquette 




1.0 


0.0 








6.89791329(1) 


Symanzik 




1.6666666 


-0.083333 








3.1814562840(7) 


TILW, p = 


S.60 


2.3168064 


-0.151791 


-0 


0128098 


1.8452250005(2) 


TILW, p = 


S.45 


2.3460240 


-0.154846 


-0 


0134070 


1.8054229585(4) 


TILW, p = 


S.30 


2.3869776 


-0.159128 


-0 


0142442 


1.7516351593(8) 


TILW, p = 


S.20 


2.4127840 


-0.161827 


-0 


0147710 


1.7188880608(5) 


TILW, p = 


S.10 


2.4465400 


-0.165353 


-0 


0154645 


1.6773505020(9) 


TILW, p = 


S.00 


2.4891712 


-0.169805 


-0 


0163414 


1.626880218(1) 


Iwasaki 




3.648 


-0.331 








0.752432061(7) 


DBW2 




12.2688 


-1.4086 








0.04881939(4) 



10 



TABLE III. Evaluation of e f 4 (cf. Eqs.(17, 19, 20), Figure 3), with Wilson gluons 
and clover fermions, for various values of the bare fermion mass m. 



m 


e4,o x 10 


e4 5 i x 10 


et± y 2 x 10 


-1.0149250 


-4.6273(2) 


1.28551(1) 


-1.85010(5) 


-0.9512196 


-4.3888(2) 


1.17807(9) 


-1.83818(4) 


-0.8749999 


-4.1089(2) 


1.05421(9) 


-1.82249(3) 


-0.8253968 


-3.9299(2) 


0.97650(7) 


-1.81156(1) 


-0.7948719 


-3.8210(2) 


0.92984(5) 


-1.80457(1) 


-0.5181059 


-2.8759(3) 


0.54981(2) 


-1.73390(5) 


-0.4234620 


-2.57089(6) 


0.43874(5) 


-1.70723(3) 


-0.4157708 


-2.54658(2) 


0.43016(6) 


-1.70500(2) 


-0.4028777 


-2.5061(2) 


0.41594(7) 


-1.70128(3) 


-0.3140433 


-2.2325(4) 


0.32325(7) 


-1.67496(2) 


-0.3099631 


-2.2202(4) 


0.31925(9) 


-1.67373(1) 


-0.3017750 


-2.1956(4) 


0.31127(7) 


-1.67125(1) 


-0.2962964 


-2.1793(3) 


0.30597(2) 


-1.66958(1) 


-0.2852897 


-2.1460(4) 


0.2953(2) 


-1.66621(2) 


-0.2825278 


-2.1379(4) 


0.2927(2) 


-1.66536(2) 


-0.2769916 


-2.1215(5) 


0.2876(1) 


-1.66366(2) 


-0.2686568 


-2.0972(1) 


0.2798(2) 


-1.66108(3) 


-0.1482168 


-1.7519(2) 


0.17636(9) 


-1.62264(9) 


0.0000 


-1.36897(4) 


0.077477(3) 


-1.57092(3) 


0.0050 


-1.35710(1) 


0.074754(4) 


-1.56906(3) 


0.0100 


-1.34534(4) 


0.072076(3) 


-1.56720(3) 


0.0140 


-1.33599(4) 


0.069967(1) 


-1.56570(2) 


0.0160 


-1.33134(4) 


0.068929(3) 


-1.56494(2) 


0.0180 


-1.32671(4) 


0.067895(2) 


-1.56419(1) 


0.0236 


-1.31382(4) 


0.065038(5) 


-1.56207(1) 


0.0270 


-1.30603(5) 


0.063335(4) 


-1.56077(3) 


0.0350 


-1.28801(6) 


0.09418(3) 


-1.55773(6) 


0.0366 


-1.28443(5) 


0.058649(4) 


-1.55711(6) 


0.0380 


-1.28131(5) 


0.057982(4) 


-1.55658(6) 


0.0427 


-1.27089(5) 


0.055767(4) 


-1.55476(7) 


0.0460 


-1.26363(5) 


0.054238(5) 


-1.55347(7) 


0.0535 


-1.24729(4) 


0.05084(2) 


-1.55051(5) 


0.0550 


-1.24405(5) 


0.05018(2) 


-1.54991(5) 


0.0720 


-1.20807(2) 


0.04295(2) 


-1.54317(6) 


0.0927 


-1.1658(2) 


0.03486(3) 


-1.53475(5) 
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TABLE IV. Evaluation of e 4 ,o (cf. Eqs.(17, 19, 20), Fig- 
ure 3), with Wilson gluons and overlap fermions, for various 
values of M , < M < 2. 



M 


e 4 ,o x 10 s 


0.01 


0.59855(5) 


0.05 


0.63347(4) 


0.10 


0.6769(2) 


0.20 


0.7628(2) 


0.30 


0.8451(2) 


0.40 


0.92220(3) 


0.50 


0.99357(4) 


0.60 


1.05872(3) 


0.70 


1.11725(2) 


0.80 


1.16893(3) 


0.90 


1.213650(1) 


1.00 


1.251396(7) 


1.10 


1.282246(6) 


1.20 


1.306376(2) 


1.30 


1.32406(2) 


1.40 


1.33561(4) 


1.50 


1.34149(2) 


1.60 


1.34201(7) 


1.70 


1.3373(2) 


1.80 


1.3271(4) 


1.90 


1.308915(8) 


1.95 


1.294780(1) 


1.99 


1.280580(7) 
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